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Derivation of the Euler Equations 
from Many-body Quantum Mechanics 



Bruno Nachtergaele^ Horng-Tzer Yau''' 

{Dedicated to Andre Verheure on the occasion of his sixtieth birthday) 

Abstract 

The Heisenberg dynamics of the energy, momentum, and particle densi- 
ties for fermions with short-range pair interactions is shown to converge to 
the compressible Euler equations in the hydrodynamic limit. The pressure 
function is given by the standard formula from quantum statistical mechan- 
ics with the two-body potential under consideration. Our derivation is based 
on a quantum version of the entropy method and a suitable quantum virial 
theorem. No intermediate description, such as a Boltzmann equation or semi- 
classical approximation, is used in our proof. We require some technical condi- 
tions on the dynamics, which can be considered as interesting open problems 
in their own right. 
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1. Introduction 

The fundamental laws of non-relativistic microscopic physics are Newton's and 
the Schrddinger equation in the classical and the quantum case respectively. These 
equations are impossible to solve for large systems and macroscopic dynamics is 
therefore modeled by phenomenological equations such as the Euler or the Navier- 
Stokes equations. Although the latter were derived centuries ago from continuum 
considerations, they are in principle consequences of the microscopic physical laws 
and should be viewed as secondary equations. It was first observed by Morrey |S] in 
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the fifties that the Euler equations become 'exact' in the Euler hmit, provided that 
the solutions to the Newton's equation are 'locally' in equilibrium. Morrey's original 
work was far from rigorous and, in particular, the meaning of 'local equilibrium' was 
not clear. It is nevertheless a very original work which led to the later development 
of the hydrodynamical limits of interacting particle systems. In terms of a rigorous 
proof along the lines of Morrey's original argument, however, there has been little 
progress until the recent work ^I]. This long delay is mostly due to a serious lack 
of tools for analyzing many-body dynamics, in the classical case and even more so 
in the quantum case. 

In this lecture, we will discuss the derivation of the Euler equations from micro- 
scopic quantum dynamics. As we want to consider the genuine quantum dynamics 
for a system with short-range pair interactions, we cannot take a semiclassical limit. 
Although one-particle quantum dynamics converges to Newtonian dynamics in the 
limit of infinite mass, this is not the case in the thermodynamic limit, i.e., the 
heavy-particle limit does not commute with the infinite- number-of-particles limit. 
This is most clearly seen in the pressure function, for which quantum corrections 
survive at the macroscopic scale. In fact, one of the conclusions of our work is that 
under rather general conditions, the pressure function is the only place where the 
quantum nature of the underlying system, in particular the particle statistics, sur- 
vives in the Euler limit. At the same time, our derivations also shows that it is the 
quantum mechanical pressure, without modification, which governs the macroscopic 
dynamics. A similar result should hold for all systems of macroscopic conservation 
laws. 

The Euler equations have traditionally been derived from the Boltzmann equa- 
tion both in the classical case and in the quantum case, see Kadanoff and Baym 0] 
for the quantum case. Since the Boltzmann equation is valid only in very low den- 
sity regions, these derivations are not satisfactory, especially in the quantum case 
where the relationship between the quantum dynamics and the Boltzmann equation 
is not entirely clear. There were, however, two approaches based directly on quan- 
tum dynamics. The first was due to Born and Green who used an early version 
of what was later called the BBGKY hierarchy, together with moment methods and 
some truncation assumptions. A bit later, Irving and Zwanzig [J used the Wigner 
equation, moment methods and truncations to accomplish a similar result. These 
two approaches rely essentially on the moment method with the Boltzmann equa- 
tion replaced by the Schrodinger equation. Unlike in the Boltzmann case, where one 
can do asymptotic analysis to justify this approach, it seems unlikely that this can 
be done for the Schrodinger dynamics. Therefore, in the present work, we follow a 
much more direct route to connect the micro- and macroscopic dynamics. 



2. Schrodinger and Euler dynamics 

We begin by considering N particles on R"^, evolving according to the Schro- 
dinger equation 

idtil^t{xi, ■ ■ ■ ,Xn) = H'iptixi, ■ ■ ■ ,xn) 
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where the Hamiltonian is given by 

^= E W{x,-x,). (2.1) 

j=l l<i<j<N 

Here, is a two-body short-ranged super-stable isotropic pair interaction and 

ipt{xi,- ■ ■ ,xn) is the wave function of particles at time t. We consider spinless 
Fermions and thus the state space is the subspace of antisymmetric functions 
in L^{M.^^). It is convenient not to fix the total number of particles and to use 
the second quantization terminology. In fact, it would be extremely cumbersome to 
work through all arguments without the second quantization formalism. The state 
space of the particles, called the Fermion Fock space, is thus the direct sum of H^: 

n := n^. 

Define the annihilation and creation operators and by 



, a; j , • • • , X]\f), 

where ax and a+ are operator- valued distributions and, as usual, " means "omit" . 

The annihilation operator ax is the adjoint of a+ with respect to the standard 
inner product of the Fock space with Lebesgue measure dx, and 

[ax, fly ]+ := axfly + axa'^ = S{x - y) , 

where 5 is the delta distribution. The derivatives of these distributions with re- 
spect to the parameter x are denoted by Va^ and Va+. With this notation, we 
can express the Hamiltonian as H = -\- V where Hq = ^ J X/a^'Vax dx and 
V = ^ J J dxdyW{x — y)ata'^ayax- It is more convenient to put the Schrodinger 
equation into the operator form, which is sometimes called the Schrodinger-Liouville 
equation. Denote the density matrix of the state at time t by 74. Only normal 
states, which can be represented by density matrices, will be considered in the 
time evolution. Then the Schrodinger equation is equivalent to idi^i = ShJi, with 
Snji '■= [H,^t] ■ The conserved quantities of the dynamics are the number of par- 
ticles , the three components of the momentum and the energy. The local densities 
of these quantities are denoted by u = {u^),lJ, = 0, • • • ,4, and are given by the 
following expressions: 

i = 1,2,3, (2.2) 
y)a+a+ayax. 

The finite volumes, denoted by A, will always by three-dimensional tori and, unless 
otherwise stated, unbounded observables on A will be defined with periodic bound- 
ary conditions. E.g., the number of particles in A, the total momentum, and the 



{ax^f{x,,---,XN) = VN + T-^''+\x,xi,--- ,xn), 
{a+^f{xu---,XM) = ^f^(-l)^-M(a.-a;,)*^-i(a;i,... 



1 If 
ut = hx = TrVa+Va^ + - dyW{x - 
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total energy of the particles in A, respectively, are defined by 



Na= dxn^, Pi^ dxp>^, j = 1,2,3, Ha = dxh^. 

J A Ja J a 

Wc slightly generalize the definition of the grand canonical Gibbs states to 
include a parameter for the total momentum of the system: the Lagrange multiplier 
a e M^. We will work under the assumption that the temperature and chemical 
potential arc in the one-phase region of the phase diagram of the system under 
consideration, such that the thermodynamic limit is unique. The infinite volume 
Gibbs states are then given by the following formula: 

r£j.XQ-l3{Ho,A + VA-Ol-PA-HN^) 

W/3,a,^(X) = ^lim^ -l^g-/3(Ho,A+VA-aPA-,.iVA) " (^"^^ 

It is convenient to denote the parameters a, /u) by A = (A'*), ;U = 0, • • • ,4 with 
A'^ = /3/U, A-' = (ia^ ,X'^=p. Define (notice the sign convention) 

3 

A-ix = Va^u" - A^w^ and {X,u)a = \M~'^ dxX{x)-u{x). 

These notations allow us to give a compact formula for the unique, translation 
invariant Gibbs state (defined with constant A) , as well as for the states describing 
local equilibrium (defined with a;-dependent A): 

u;x= lim el^l<^''*>^/ZA(A) (2.4) 

where Za{X.) is the partition function given byZA(A) = Trel^l<^'">A. If we define the 
pressure, as a function of the constant vector A, by il>{X) = limi,_>oo | A|~^ log Za{X), 
then 

^ = -a(u-). (2.5) 

As the states u)\ are translation invariant, we have u)x{u'^) = q^. Explicitly, 

P = WA(nx), c\=uj\{px), e = ux{hx). 

Notice that q and e are momentum and energy per volume. 

Again, we will work under the assumption that these parameters stay in the 
one-phase region, the limiting Gibbs state is unique and these definitions are un- 
ambiguous. Although momentum is preferable as a quantum observable, we also 
introduce the velocity in order to be able to compare with the classical case. The 
velocity field v{x) has to be defined as a mean velocity of the particles in a neigh- 
borhood of X. Therefore we have v{x) = q{x)/p{x). We also introduce the energy 
per particle defined by e = e/p. The usual Euler equations are written in terms of 
p, V, and e. 

In order to derive the Euler equations, we need to perform a rescaling. So we 
shall put all particles in a torus A^ of size and use {X,T) = {ex,et) to denote 
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the macroscopic coordinates. For all equations in this paper periodic boundary 
conditions are implicitly understood. 
The Euler equations are given by 




These equations are in form identical to the classical ones but all physical quantities 
are computed quantum mechanically. In particular, P{e,p) is the thermodynamic 
pressure computed from quantum statistical mechanics for the microscopic system. 
It is a function of X and T only through its dependence on e and p. If no velocity 
dependent forces act between the molecules of the fluid under consideration (we 
consider only a pair potential) , the pressure is independent of the velocity. 
Let q = {q'^, • • • , q*), related to density, momenta and energy as follows: 

q" = p, q'=pv\ q^^e^pe. (2.7) 

In other words q^,q^,q^ , and q'^ arc momenta and energy per volume instead of per 
particle as in the usual Euler equation H2.6fl . We rewrite the Euler equations in the 
following form 

^ + ^Vf [4^9)] =0, fc = 0,l,2,3,4. (2.8) 

i—l 

The matrix A is determined by comparison with the Euler equations: 
A° = q^, A) = S^jP + q^qJ/qo, A] = q\q^ + P)/go. 



3. Local equilibrium 

To proceed we need a microscopic description of local equilibrium. Suppose 
we are given macroscopic functions qiX). We wish to find a local Gibbs state with 
the conserved quantities given by q{X). The local Gibbs states are states locally in 
equilibrium. In other words, in a microscopic neighborhood of any point x G the 
state is given by a Gibbs state. More precisely, we wish to find a local Gibbs state 
with the expected values of the energy, momentum, and particle number per unit 
volume at X given by q{X). To achieve this, we only have to adjust the parameter 
A at every point X. More precisely, we choose X{X) such that the equation (|2.5|l 
holds at every point, i.e.. 



5*(A(X)) _ 



^X^'{X) 



q^{X). 
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If we denote the solution to the Euler equation by q{X,T), then we can choose in 
a similar way a local Gibbs state with given conserved quantities at the time T. 
Define the local Gibbs state 

^1 = ^ exp [e-^(Xiet, e ■ ), (3.1) 

where Ce(t) is the normalization constant. Clearly, we have that uji;{t){ui^) = 
q'^{ex,et) to leading order in e. 

In summary, the goal is to show that, in the limit e ^ 0, the following diagram 
commutes: 



q(X,0) q{X,T) 



local equilibrium 



limit e I of expectation 
of locally averaged observ- 
ables 



Schrodinger 

70 > 



As smooth solutions of the Euler equations are guaranteed to exist only up 
to a finite time [S], say Tq, we will formulate our assumptions on the dynamics of 
the microscopic system for a finite time interval as well, say t G [0,To/e]. Note the 
cutoff assumptions below would hold automatically for lattice models. 



4. Assumptions and the main theorem 

Our main result is stated in Theorem 14. II below. First, we state the assump- 
tions of the theorem with some brief comments. There are three kinds of assump- 
tions. 

The first category of assumptions could be called physical assumptions on the 
solution of the Euler equations that we would like to obtain as a scaling limit of the 
underlying dynamics, and on the pair interaction potential of this system. 

I. One-phase regime: We assume that the pair potential, W, is and with 
support contained in a ball of radius R. Moreover, we assume that W is symmetric 
under reflections of each of the coordinate axes (e.g., rotation symmetric potentials 
automatically have this symmetry), and has the usual stability property |12| : there 
is a constant B > such that, for all > 2, xi, . . . , G M.^, 

W{xi-Xj)>-BN. 

l<l<j<N 

Of the Fermion system with potential W we assume that there is an open 
region _D C K.^, which we will call the one-phase region, such that the system has 
a unique limiting Gibbs state for all values of particle density and energy density 
ip,e)GD. 

The solution of the Euler equations we consider, q{X,T), will be assumed to 
in A" for r e [0, To], and have local particle and energy density in the one-phase 
region for all times T e [0,ro]- I-c, ipiX,T), e{X,T)) e D, for all X € Ai and 

re [o,To]. 
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The next category of assumptions is on the local equilibrium states for the 
Fermion system that we construct and on their time-evolution under the Schrodinger 
equation. They can be considered conjectures. In fact, these assumptions have not 
been rigorously proved even for Gibbs states. Although one expects that these 
assumptions can be proved in the high-temperature and low-density region using 
some type of cluster expansion methods, this has only been done recently for Bosons 
in [7]. For the rest of this paper, we shall assume this cutoff assumptions for the 
solution to the Schrodinger equation that we consider, as well as for the Gibbs states 
in the one phase regions considered in this paper. 

II. Cutoff assumptions: Suppose that 7t is the solution to the Schrodinger equa- 
tion with a local equilibrium state as initial condition, constructed with the param- 
eters derived from a solution of the Euler equations (with the appropriate pressure 
function) that does not leave the one-phase region. We make the following two 
assumptions about this solution: 

1. Finite velocity cutoff assumption: Let Np{t) ~ Tr7ta+ap, where is the Fourier 
transform of af. Then there is a constant c > 0, such that for all t < Tq/e, 



2. Non-implosion assumption: There is a constant Ctq, such that for all t < Tq/e, 



where R is the range of the interaction W. 

Finally, we have an assumption on the set of the time-invariant ergodic states 
of the Fermion system. To state this assumption we need the notion of relative 
entropy, of a normal state 7 with respect to another normal state lu. Let 7 and 
ui denote the density matrices of these states. The relative entropy, S'(7 | w), is 
defined by 



For a pair of translation invariant locally normal states, one can show existence of 
the relative entropy density JOIi defined by the limit 



where 7a^ and uj\^ denote the density matrices of the normal states obtained by 
restricting 7 and to to the observables localized in = £~^Ai. The existence of 
the limit can be proved under more general conditions on the finite volumes, but 
this is unimportant for us. 

III. Ergodicity assumption ("Boltzmann Hypothesis"): All translation in- 
variant ergodic stationary states to the Schrodinger equation are Gibbs states if 
they satisfy the following assumptions: 1 ) the density and energy is in one phase 
region. 2) The relative entropy density with respect to some Gibbs state is finite. 





' 



(4.1) 




if kcro; C ker7, 
else. 



s(7|w) = lim£^5(7A^ | waJ , 
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We expect that our assumptions hold for the solutions 7t of the Schrodinger 
equation that we employ, but it must be said that, at this moment, very little is 
known. We believe that these are natural conditions. To prove them under rather 
general conditions or even for a special class of models may be regarded as an 
important open problem in quantum statistical mechanics. 

Theorem 4.1 Suppose that q{X,T) is a smooth solution to the Euler equation in 
one phase region up to time T <Tq. Let ujf be the local Gihhs state with conserved 
quantities given by q{X,T). Suppose that the cutoff assumptions and the ergodicity 
assumption hold. Let 74 be the solution to the Schrodinger equation and 70 — lOq 
(Notice 7t depends on e). Then for all t < To/e we have 



In other words, cuf is a solution to the Schrodinger equation in entropy sense. In 
particular, for any smooth function f on A, we have 



This theorem is a quantum analogue of the classical result by OUa, Varadhan, 
and Yau There are a few differences in the assumptions and the strategy 

followed in their paper with respect to ours. E.g., in the treatment of of 
the classical case the cut-off assumption 1) was not needed. Instead, the usual 
quadratic kinetic energy was replaced by one with bounded derivatives with respect 
to momentum. For Fermion models on a lattice instead of in the continuum, no cut- 
off assumptions are required. Another difference with the treatment in is that 
we do not add noise terms to the "native" dynamics. In [TT) a weak noise term was 
added to the Newtonian dynamics in order to be able to prove convergence to local 
equilibrium. In that paper, the strength of the noise vanishes in the hydrodynamic 
limit. Here, we do not modify the Heisenberg dynamics in any way, but instead 
reduce the question of convergence to local equilibrium to the ergodicity property 
given in Assumption III. 

5. Outline of the proof 

The basic structure of our proof follows the relative entropy approach of [111 
I13j . The aim is to derive a differential inequality for the relative entropy between 
the solution to the Schrodinger equation and a time-dependent local Gibbs state 
constructed to reproduce the solution of the Euler equations. The time derivative 
of the relative entropy can be expressed as an expectation of the local currents with 
respect to the solution to the Schrodinger equation: Since 74 is a solution to the 
Schrodinger equation, we have for any density matrix lu^ the identity 



lim s(7t|wf) = 0. 




^SijtH) = Tr7t {i6h ~ dt}\ogLot. 



(5.1) 
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Using the definition of the local equilibrium states H2.4|l . we get 

^s(7tkD = £'Tr7t {z5h - 9J {(A,(i, •), u) - logc,(i)} • (5.2) 

Direct computation yields that (5// (A£(t, •), u) — — eX]j=i(^j'^e(^i ')i ^j(^))' to 
leading order in e, and where the 6j are local observables for the microscopic 
currents of the conserved quantities. Since we do not know the solution of the 
Schrodinger equation, 74, well enough, this expectation in the RHS of (|5.2|l cannot 
be computed or estimated explicitly. The main idea is to bound these expectations 
in terms of the relative entropy itself, thus obtaining a differential inequality. All 
these bounds are based on the following inequality, which is a consequence of the 
variational principle for the relative entropy [HI llOj : for density matrices/states to 
and 7, with kerw — {0}, and for all self-adjoint /i, and (5 > 0, one has 

-f{h) < 6-^ logTr e*''+i°g- + d-^S{-f\Lu). (5.3) 

Most of the work is then to show that the terms resulting from the first term in the 
RHS of this inequality are sufficiently small. This requires a number of steps. Two 
essential ingredients are the Euler equations (naturally), and a quantum version 
of Virial Theorem to relate certain terms in the currents to the thermodynamic 
pressure as given by quantum statistical mechanics. Next, we explain the main 
steps in a bit more detail. 

Step 1: Replace the local microscopic currents by macroscopic currents. The 
basic idea in hydrodynamical limit is first to show that the local space time average 
of the solution is time invariant. From Assumption III, ergodic time invariant 
states are Gibbs. For Gibbs states, we can replace the local microscopic currents 
by macroscopic currents. 

la: Construct a commuting version of local conserved quantities. In principle, 
one should be able to express the macroscopic currents as functions of the con- 
served quantities. The local conservative quantities, however, are operators which 
commute only up to boundary terms. Therefore, we either need to prove that the 
non-commutativity does not affect the meaning of macroscopic currents or we need 
to construct some commuting version of the local conservative quantities. We fol- 
low the second approach and construct a commuting version of local conservative 
quantities with a method inspired by 

lb: Restriction to the one phase region. Since Assumption III can only be 
applied in the one phase region, we have to introduce suitable cutoff functions. As 
these, in general, do not commute with the local currents, this is a non-trivial step. 

Ic: Apply the Virial Theorem. Following the method of [Sj, we prove and 
apply a quantum virial theorem to relate the local currents to the pressure. 

Step 2: Estimate all errors by local conservative quantities. For this step the 
entropy inequality (|5.3|l is crucial. 

Step 3: Derive a differential inequality of the entropy with error term given by 
a large deviation formula. As there is no large deviation theory for non-commuting 
observables, it is essential that we have expressed everything by commuting objects. 
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The rest of the argument follows by the standard relative entropy method. 
The main technical difficulties, in comparison with the classical case, all stem from 
the non-commutativity of the algebra of observables. Simple inequalities, such as 
|A + -B| < |^| + |-B| and \AB\ < \A\\B\, which are used numerous times in estimates 
for classical systems, are false for quantum observables. Therefore, all estimates 
have to be derived without taking absolute values. A full account of our work will 
appear shortly |T1] . 
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